Statistical properties of the critical eigenstates in power-law random banded matrices 

across the band 
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The level-spacing distribution in the tails of the eigenvalue bands of the power-law random banded 
matrix (PRBM) ensemble have been investigated numerically. The change of level-spacing statistics 
across the band is examined for different coupling strengths and compared to the density of states for 
the different systems. It is confirmed that, by varying the eigenvalue region, the same level-spacing 
statistics can be reached as by varying the coupling strength. 
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The Anderson metal-insulator transition (MIT) is a 
phenomenon of major physical importance that contin- 
ues to attract a substantial research effort^. With just 
short-range interactions, localization-delocalization (LD) 
transitions are only found in systems with dimension- 
ality, D, greater than two4 However, with the addi- 
tion of long-range interactions, or correlations between 
the short-range interactions, it is possible to study the 
LD transition in systems with dimensionality less than 
two*. In this respect, power-law random-banded matri- 
ces (PRBMs), that exhibit this transition, have recently 
attracted much attention 6 : 7 : 8 : 9 : 10 . The PRBM ensemble 
was introduced by Mirlin et. aliii and, in the real case, 
is defined as the ensemble of N x N random symmetric 
matrices, H. The PRBM elements, Hij, are randomly 
drawn from a Gaussian distribution, centred around zero, 
with a variance governed by a power-law decay: 
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where a and b € (0, oo) are parameters. Regarding H as 
a Hamiltonian, the eigenvalues, E, are energies. 

For a = 1, it has been shownii that all the eigen- 
states of these matrices are critical (i.e. at the LD tran- 
sition). The parameter b is inversely related to the cou- 
pling strength between the nodes. In the limit b 3> 1 and 
a = 1, the PRBM critical states are analogous to the 
critical states at the Anderson transition with D = 2 + e 
and e <C 1, and for b <C 1 to those found in the Ander- 
son model with D ^> 1—. By varying b, it is possible 
to access a set of different critical theories parameter- 
ized by dimension (2 < D < oo) in the conventional 
Anderson transition. This ability to examine Anderson 
transitions in different dimensions, in the same effectively 
one-dimensional model, makes the study of PRBMs a 
powerful method to make progress in this rich field. As 
well as being an analogue for the study of important tran- 
sitions elsewhere, the PRBM is physically important in 
its own right and has been applied to the study of the 
finite-temperature Luttinger liquid^, the coherent prop- 
agation of two interacting particles in a ID weak ran- 
dom potential^ and other problems 7 ^. Recently, it has 



also been realized that, with the addition of chirality, the 
PRBM also describes the LD transition of quark zero 
modes in QCD 8 ". In the QCD vacuum, the quark zero- 
mode wavefunction decay has a power-law dependence 
and long-range hopping between sites is possible. In this 
model, the eigenvalues away from the centre of the spec- 
tral band are not affected by the chiral structure; this 
makes it relevant to the analysis of the critical states 
away from the band centre given in this paper. 

It is known— that, by changing the coupling parame- 
ter, b, for a — 1, it is possible to access the set of mod- 
els in different dimensions by investigating the statistical 
properties of the PRBM eigenstates, E, around the band 
centre, E = 0. Alternatively, this can be achieved by 
varying E away from the band centre and keeping b con- 
stant, as was mentioned in RefA However, this feature 
of the PRBM has not been previously studied and con- 
firmed numerically, and this is the aim of this paper: to 
investigate the statistical properties of the eigenstates of 
the PRBM across the whole range of the spectrum. One 
of the main results of the paper is to demonstrate numer- 
ically that there is a mapping between different values of 
6 and E given certain critical level-spacing statistics, i.e. 
(\b'\,\E'\) ^ {\b"\,\E"\) with \b"\ > \b'\ and \E"\ > \E'\, 
in terms of equivalent level-spacing statistics. This is im- 
portant, both in models where PRBMs are used to model 
systems in which all the eigenvalues are of physical im- 
portance, see e.g. Ref. 8 , and to support the validity of 
taking statistics from a finite region iiSiiS around the 
band centre; we obtain estimates to justify the width of 
the region around the band centre that is used for statis- 
tical analysis. 

To study the level-spacing statistics, a careful unfold- 
ing procedure was used, taking into account the varying 
mean level spacing (especially in the band tails) when 
normalizing the spacings. The width of the bins in E 
was chosen to be sufficiently small so as to ensure that the 
nature of the states, and thus their statistical properties, 
do not change appreciably across them. The symmetry 
of the spectrum about E = has also been used to im- 
prove the level-spacing statistics. The size of the N x N 
matrices (with periodic boundary conditions) varied in 
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FIG. 1: Density of states of the PRBM ensemble with a = 1 
and varying values of b as marked. The inset shows portions of 
the densities of states after scaling so that 90% of the states 
are within the interval E' G [—1,1]. The arrows show the 
eigenvalue position corresponding to the 95th percentile of 
the DOS. 

the range between N = 1000 and N = 8000. The off- 
diagonal elements of the matrices were multiplied by a 
factor 0.5 (as in Ref. — ) in order to reinforce the close- 
ness of the system to the Gaussian orthogonal ensemble. 

The densities of states (DOS) of the PRBM ensembles 
for different values of the parameter b and a = 1 are 
shown in Fig. ^ It can be seen that with increasing b the 
bandwidth increases. There is substantial deviation from 
the semi-circle shape (cf. Ref. — ) typical for the spec- 
tra of matrices from the Gaussian orthogonal ensemble, 
which do not show any change in the level-spacing statis- 
tics (the Wigner surmise) across the band— The inset in 
Fig-Hshows eigenvalue-rescaled DOSs, and verifies what 
is suggested by the main graph - that lower values of b 
also correspond to longer tails appearing in the density 
of states. This means that the DOSs corresponding to 
different values of b cannot be equated by a simple linear 
rescaling of the eigenvalues. The fact that the plots in 
the inset appear to intersect at one point is not signifi- 
cant, but merely a result of the particular variables used 
to produce the graph. 

The bands shown in Fig. ^ contain only critical states 
of multifractal nature, which exhibit a particular type 
of level-spacing statistics. For these critical states, the 
probability distribution function, P c (s), of normalized 
level spacings, s, between adjacent levels starts from zero 
at s — (level-repelling effect), goes through a maximum, 
and then exponentially decays for s>l according to the 
following law&i^: 

P c (s)oce- As " . (2) 

This functional form is different from those appropri- 
ate for descriptions of level-spacing statistics for local- 



FIG. 2: Double logarithmic plots of — ln/(s) versus s for 
a = 1, 6 = 1, in the central eigenvalue region \E\ G [0,0.3] 
with a straight-line fit (solid line) to the displayed region, 
giving aslope f3 = 1.58 and the band-tail region \E\ G [1.6, 1.7] 
(dashed line), yielding a gradient of 1.38. Configurational 
averaging over 10 5 realizations of 1000 x 1000 matrices was 
used for the statistical analysis. The inset shows how f3 varies 
with b at the band centre obtained by linear fits in the range 
s G [1.5, 3] (crosses) with the similar dependence obtained 6 
by linear fits for s > 2 (circles). (The curves in the inset are 
guides to the eye.) 



ized or extended states. The level spacing for local- 
ized states follows the Poisson law, Pp(s) = exp (— s), 
while for extended states, the distribution is described 
by the Wigner-Dyson surmise distribution—, Pw(s) = 
(7r/2)sexp(-(7r/4)s 2 ). 

In the case of the standard Anderson transition for 
systems in Euclidean geometries, there has been some 
controversy about the value of the exponent f3 in Eq. J2J. 
Recent papers&21 have given (3 ~ 1.2 for 3D, whilst this 
value was refuted in another paper^ and the value set at 
(3 = 1.0 ± 0.1. The uncertainty in the numerical value of 
f3 for the standard Anderson transition comes from the 
fact that it is difficult to evaluate this exponent from an 
analysis of the critical states which exist only around one 
point (the LD transition) on the energy scale. This dif- 
ficulty is lessened for the PRBM where all the states are 
critical, and thus an analysis of the level-spacing statistics 
can be performed relatively easily under the assumption 
that the statistics do not change with varying E. This 
assumption has been made 6 to perform a level-spacing 
statistical analysis around E = Q. It was demonstrated 
there that the exponent (3 in Eq. (0) depends on the pa- 
rameter b. An obvious question, which has not been ad- 
dressed before, is how the level-spacing statistics change 
away from the band centre. In particular, is there a cor- 
respondence between the statistics for different values of 
E and 6? These questions are answered below. 

In order to examine the level-spacing statistics for dif- 
ferent values of E, we have performed a standard anal- 
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FIG. 3: The exponent (3 versus the eigenvalues of the PRBM 
for different values of parameter b (marked by different sym- 
bols). The arrows indicate the eigenvalues corresponding to 
the 95th percentile of the respective DOS (see Fig. 0. Two 
data sets for b — 1 (open diamonds obtained for 1000 x 1000 
matrices averaged over 10 5 configurations, and solid diamonds 
obtained for 8000 x 8000 matrices averaged over 420 configu- 
rations) show the absence of a size effect. The bin width on 
the eigenvalue axis was independent of eigenvalue and equal 
to 0.1 for all the models. The symbols, * and #, denote two 
pairs of parameters (\E\, b) for which the level-spacing statis- 
tics have the same exponent (3. 



ysis, studying the cumulative level-spacing distribution 
I(s) = f P(s')ds', which retains the asymptotic be- 
haviour of P(s). We expect a linear dependence in a 
double- logarithmic plot of ln(J) vs s for s ^> 1, with the 
gradient (i.e. (3) being in the range (3p < (3 < /?Wi where 
f3p = 1 and /3w = 2 are the exponents for the Poisson 
and Wigner-Dyson distributions, respectively. Indeed, 
we have found such linear dependences characterized by 
different values of (3 for varying E. In Fig.0 we show two 
representative dependencies for the states at the band 
centre (solid line, (3 = 1.58) and from the tail region 
(dashed line, (3 = 1.38). We should emphasize that we 
have used a slightly different statistic from that of Ref£. 
This was due to the fact that the DOS in the tail regions 
rapidly approaches zero and there are a smaller number 
of eigenvalues (and thus level spacings) available for sta- 
tistical analysis. Therefore, we have moved the range of 
level spacing s towards smaller values, s 6 [1.5, 3], where 
there are more data in P(s) and the double-log graph 
is still reasonably straight, in contrast to the higher s 
ranges used in Ref. 6 «. Such a change in the range of anal- 
ysis resulted in slightly higher values of (3 for the states 
around the band centre (see the inset in Fig. |2J . We have 
checked our method of analysis for the band centre using 
higher s ranges and reproduced the published values^ of 
(3 for different b. 

Fig.|3|shows our main result. In this figure, we demon- 
strate the possibility of mapping between different values 



FIG. 4: The deviation, A/3 = (3(E) - /3(0), from its value at 
the band-centre, plotted for different values of b. The eigen- 
value abscissa has been scaled to give the percentile of the 
density of states, with the band centre taken as zero and the 
far tail as 100%. 



of eigenvalue E and parameter b in terms of the same 
level-spacing statistics characterized by the exponent (3. 
Indeed, as can be seen from Fig. [31 one value of \E\ = \E'\ 
for a given b = b' (e.g. the star-symbol) corresponds to a 
different value of | £7 1 = \E"\ with another value of b = b" 
(e.g. the hache symbol) but having the same statistics 
characterized by (3. In other words, we demonstrate that 
it is possible to access different dimensionalities for the 
standard Anderson LD transition simply by analyzing 
different eigenvalues of the PRBM at fixed value of the 
the parameter b. 

The plots in Fig. [3] clearly show a region around the 
band centre where the statistics do not change much 
followed by a decrease in (3 towards the Poisson value 
((3 = 1) on approaching the edge of the band. The 
range of eigenvalues around the band centre where the 
value of (3 is approximately constant depends on the 
value of parameter b. In order to determine the extent 
of this region around the band centre, it is illuminating 
to plot the change of the exponent (3 from the band cen- 
tre, A/3 = (3(E) — (3(0), versus the cumulative density 
of states, N(E) = 2 g(E)dE. In this case, as seen 
in Fig. all the curves in the midgap region, say for 
N(E) < 0.5, approximately collapse to a single curve. 
This justifies an analysis of the level-spacing statistics 
around the band centre for less than approximately 50% 
of states with an absolute error in (3 less than 0.1. 

The data presented in Fig. |31 were obtained for 8000 x 
8000 matrices by averaging over 420 realizations for all 
values of b except 6 = 2, for which 620 realizations were 
used to obtain better statistics in the far tail. The ab- 
sence of size effects has been checked for the case of b = 1 
by performing additional analyses for 1000 x 1000 matri- 
ces (compare open and solid squares in Fig. [3Jl . 



4 




1 2 3 4 5 



.1 



FIG. 5: The n-level spacing probability distribution func- 
tions P(n,s) (n = 0, ...,4) for two sets of parameters (b 1 = 
0.6, \E'\ e [0,0.3]) and (6" = 1, \E"\ G [1.6,1.7]), both char- 
acterized by the same value of j3 ~ 1.41. For generating these 
graphs, 1000 x 1000 matrices averaged over 25000 (&' = 0.6) 
and 10 5 (b" = 1) realizations were used. 

In order to verify the equivalence between different val- 
ues of (6, \E\) in terms of equivalent level-spacing statis- 
tics, found by analysis to have the same value of the expo- 
nent (5 for the adjacent level-spacing probability distribu- 
tion P c (s), given by Eq. (2), we have calculated in addi- 



tion the n-level distribution functions, P c (n, s)i&. These 
functions are the probability distribution functions for 
the level spacing s between two levels separated by n 
other levels, with P c (n = 0, s) = P c {s)- Fig. shows the 
equivalence between all five n-level distribution functions 
(n = 0, . . . , 4) for two representative pairs of parameters, 
(6', l-E'l) and (6", \E"\) characterized by the same value 
of 1.41). This means that the statistical proper- 
ties of the critical eigenvalues, which can be character- 
ized by different functions, e.g. the two-point correlation 
function—, B, c (s) = S(s) — 1 + J2nLi Pc(n, s), are the 
same for (b',\E'\) and (V, \E"\). 

In conclusion, we have provided numerical evidence for 
the equivalence of statistical properties of critical eigen- 
states for different values of both the strength parameter 
b and eigenvalues of power-law random banded matrices 
(PRBMs). This justifies the use of an alternative way 
of exploring different dimensionalities for the standard 
Anderson transition. Instead of the usual way of ana- 
lyzing the level-spacing statistics around the band centre 
for different values of b, we show that it is possible to 
access the same statistics by moving through the band 
of eigenvalues at fixed value of b. Also, we specify the 
range of the eigenvalues around the band centre where 
the level-spacing statistics do not change significantly, 
and thus justify the standard method of statistical anal- 
ysis of eigenvalues of the PRBM by varying b. 
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